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Abstract 

We consider the 2d XY Model with topological lattice actions, which are 
invariant against small deformations of the field configuration. These ac- 
tions constrain the angle between neighbouring spins by an upper bound, 
or they explicitly suppress vortices (and anti- vortices). Although topo- 
logical actions do not have a classical limit, they still lead to the universal 
behaviour of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition 
— at least up to moderate vortex suppression. Thus our study un- 
derscores the robustness of universality, which persists even when basic 
principles of classical physics are violated. In the massive phase, the an- 
alytically known Step Scahng Function (SSF) is reproduced in numerical 
simulations. In the massless phase, the BKT value of the critical expo- 
nent Tjc is confirmed. Hence, even though for some topological actions 
vortices cost zero energy, they still drive the standard BKT transition. In 
addition we identify a vortex-free transition point, which deviates from 
the BKT behaviour. 
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1 Introduction 

Universality is of central importance in quantum field theory and statistical me- 
chanics, because it makes the long-distance physics insensitive to the short- distance 
details at the cut-off scale. The corresponding universality classes are determined 
by the space-time dimension and by the symmetries of the relevant order parame- 
ter fields. In lattice field theory, one often demands, in addition, the lattice action 
to have the correct classical continuum limit. Recently, we have introduced the 
concept of topological lattice actions, which do not have a classical limit Topo- 
logical lattice actions are invariant against small deformations of the lattice fields. In 
0{N) Models, the simplest topological action constrains the relative angle between 
nearest-neighbour spins to a maximal angle 6. All allowed configurations (that do 
not violate this constraint) are then assigned the action value zero. Since the ac- 
tion does not vary at all, it does not give rise to a meaningful classical equation of 
motion. Consequently, it does not have the correct classical continuum limit, and 
perturbation theory does not apply either. As we have demonstrated analytically 
for the Id 0(2) and 0(3) Model, despite this classical deficiency, the topological 
lattice action still leads to the correct quantum continuum limit. However, for these 
Id topological actions the lattice artifacts go to zero only as 0{a) in the limit of 
vanishing lattice spacing a, while they are of O(a^) for the standard lattice action. 

The correct quantum continuum limit has also been verified in the 2d 0(3) 
Model Based on numerical simulations with the Wolff cluster algorithm 0], we 
have reproduced the analytic results for the Step Scaling Function (SSF) |3| that 
was introduced in Ref. j^]- Remarkably, in the well accessible range of correlation 
lengths, the cut-off effects of the topological action are smaller than those of the 
standard action and of the tree-level improved Symanzik action, which had been 
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investigated previously jsf. By combining the standard and the topological action, 
we have constructed a highly optimised constraint action for the 2d 0(3) Model that 
has only per mille level cut-off effects of the SSF for ratios a/L < 0.1 [6]. Although 
the topological susceptibility receives contributions from zero-action dislocations, it 
was found to diverge only logarithmically rather than with a power law, as a 
semi-classical argument would suggest While it has been suspected that 9 is an 
irrelevant parameter which gets renormalised non-perturbatively, we have identified 
distinct physical theories for each value < < vr [s*] (see also Refs. (9|,[l3). At 
^ = we also investigated a topological lattice action which explicitly suppresses 
topological charges. Although this action does not have the correct classical contin- 
uum limit either, it was found to have the correct quantum continuum limit as well 

Q. 

This paper addresses the 2d XY (or 0(2)) Model, which has been applied, for 
instance, to describe thin films of superfluid helium fluctuating surfaces and 



their roughening transition, as well as Josephson junction arrays [12|. Here we inves- 
tigate topological lattice actions for that model. In contrast to the 2d 0(3) Model, 
which is asymptotically free, the continuum limit of the standard 2d XY lattice 
model is reached at finite values of the coupling. It corresponds to the well-known 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition, an essential transition of 



infinite order {131 . Il4j |. The BKT transition separates a massive phase, in which 
vortices are condensed, from a massless phase, with bound vortex-anti- vortex pairs. 

Although it is not asymptotically free in the usual sense, the 2d XY Model has a 
non-trivial massive continuum limit at the BKT phase transition. There is numerical 
evidence that this continuum limit corresponds to the sine-Gordon Model at coupling 



/3 — VSvr [15|, which in turn is equivalent to the SU{2) chiral Gross-Neveu Model. 
In this sense, the continuum theory is asymptotically free after all. The SSF 
has been worked out analytically, and tested against numerical simulations [if 
Remarkably, in this case even the cut-off effects, which vanish only logarithmically 



as one approaches the continuum limit, have universal features [17 . 

It is interesting to investigate whether topological lattice actions lead to the 
usual quantum continuum limit also in this case. One question is how far universality 
really reaches, in view of the critical behaviour, and of the cut-off effects. As a further 
motivation, we refer to an estimate of the critical temperature for the standard 
lattice action, based on the energy cost for isolated vortices (or anti- vortices), which 
tend to disorder the system. If this is a relevant argument behind the BKT phase 
transition, then the behaviour for topological lattice actions is in fact tricky. 

Some time ago, the BKT phase transition has been investigated in the so-called 



Step Model [18H21| . The Step Model has a topological action, which vanishes if the 
angle between nearest-neighbour spins is less than 7r/2; otherwise it is a positive 
constant 5*0 El While in the Step Model the BKT transition is attained by varying 
5*0, it is attained with the constraint action by varying 6. As 5*0 is sent to infinity. 



"'^Also the version with a finite step at a variable angle has been addressed with analytical 
approaches 1^, 20j . 
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the Step Model approaches the constraint action with 5 = 7r/2. On a square lattice, 
vortices are completely eliminated in that case. In agreement with the BKT picture, 
this point in the phase diagram turns out to be in the massless phase. For smaller 
values of 5*0, vortices have a finite action. After some controversy, it has been 
confirmed that the Step Model is indeed in the BKT universality class 221-124 . 



Using efficient cluster algorithms, we will show in this paper that the constraint 
angle action also falls into the BKT universality class. This follows by comparison 



with analytic results for the SSF j25|, l26|, and for the critical exponent r^c j27|, |28 



On the other hand, the cut-off effects of this topological action do not share the 
predicted universal features. 

We further investigate a topological action that combines the constraint angle 6 
with explicit vortex suppression, by assigning an action value A > to each vortex 
or anti-vortex. Also that action turns out to have the universal features of the BKT 
transition, at least up to A ~ 4. A different behaviour is observed, however, at 
the endpoint of this transition line, which seems to be located at S = n (no angle 
constraint) and A = oo (no vortices). 

In Section 2 we describe topological actions with two parameters, for an angle 
constraint and an explicit vortex suppression. Section 3 investigates these actions — 
with the angle constraint included — by approaching the phase transition both in 
the massive and in the massless phase. In Section 4 we address a topological vortex 
suppression action without an angle constraint, and the extrapolation A — ?■ +oo. 
Section 5 contains our conclusion. Finally the cluster algorithm used to simulate the 
topological actions is explained in Appendix A, and Appendix B discusses surprising 
aspects of the correlations in ferromagnetic systems. 



2 Topological Lattice Actions 

Let us consider the 2d XY Model on a periodic square lattice. A 2-component unit 
vector Cx = (cosfyO^;, sinf/}^,) is attached to each lattice site x. The standard lattice 
action reads 

5'standard[e] = - 6^ ■ Cy] = ^1 - COs{Lpx - ^fy) , (2.1) 

{xy) (xy) 

where (xy) denotes a pair of nearest-neighbour sites, and the parameter (3 corre- 
sponds to an inverse coupling. A vortex number va G {0,±1} is associated with 
each elementary plaquette □, with the corners Xi, X2, X3, x^ in counter-clockwise 
order. Introducing the relative angles 

^^{x,xj) = {fx, - fxj) mod 27r G (-vr, n] , (2.2) 

the vortex number of a plaquette is given by 

= ^ {^V{xix2) + ^V{x2xz) + ^^{xzx^) + ^^{x^x^)) e {0, ±1} . (2.3) 
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Higher vortex numbers cannot occur. The vortices are known to be the relevant de- 
grees of freedom that drive the BKT phase transition 1^- According to Stokes' 
Theorem, the sum of all vortex numbers on a periodic lattice always vanishes, 



Let us now introduce a topological action as a sum over elementary plaquettes, 

S[e] = xY,\vo\ . (2.4) 



□ 



This action counts the number of vortices (with va = 1) plus anti-vortices (with 
fn = —1), and multiplies this sum with the single-vortex action A. In particular, 
the limit A — >■ oo removes all vortices. When one continuously varies the spin field, 
without changing the (discrete) vortex number \va\, the action does not change 
either. Consequently, it is invariant against small deformations of the lattice field, 
so it represents a topological action. 

Let us mention that the analogous A-term has also been introduced in the 3d XY 
Model 29| and 0(3) Model j30|. In both cases it was combined with the standard 



term to investigate the phase diagram with the axes /3 and A. This also involved 
studies of the topological action at /3 = 0, where phase transitions at finite Ac were 
observed. 

We may further modify the pure vortex suppression action by imposing the 
angle constraint |Ay9^a,j^)| < S, which restricts the relative angle Aipi^^y) between 
nearest-neighbour spins Cx and Cy to a maximal value 6 G [0, vr]. Allowed configura- 
tions (which obey this angle constraint) still have the action value S[e] of eq. fl2.4p . 
while all other configurations (which violate the constraint on at least one nearest- 
neighbour pair of sites) are assigned an infinite action, so they are eliminated. The 
actions characterised by the parameter A and the angle constraint 6 remain invariant 
under small field deformations, and are thus still topological. 



3 Universal Behaviour of Angle Constraint Topo- 
logical Actions 

In this section, we investigate the 2d XY Model with topological lattice actions that 
impose an angle constraint for nearest-neighbour spins, S < n. In addition, the 
actions may or may not explicitly suppress vortices, A > 0. The universal behaviour 
is studied both in the massive and in the massless phase. 



3.1 Phase Diagram 

To determine the critical angle 6c of the constraint topological action, we measure 
the correlation length ^(5) in the massive phase close to the phase transition that 
occurs in the infinite volume limit. This is done by increasing the lattice volume 
V = L X L until the correlation length ^{6, L) converges to its infinite volume limit. 
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For angles 5 > 5c, not too close to the phase transition, the convergence is observable 
on tractable lattice sizes (up to L = 2000). To determine the critical point 5c we 
fit the correlation length S,{5) to a function, which is characteristic for the BKT 
transition, 



5 -5c 



l/2^ 



(3.1) 



where A and B are fitting parameters. This form represents an essential {i.e. infinite 
order) phase transition (for conventional lattice actions, the coupling 1/a/5 takes 
the role of 5). The critical angles 5c obtained from these fits (which have a good ratio 
X^/d.o.f.) are listed in Table [T]for the topological action without vortex suppression, 
A = 0, and with explicit vortex suppression, A = 2 and A = 4. 



A 


5c 





1.77521(57) 


2 


1.86648(81) 


4 


1.9361(83) 



Table 1: Critical angles 5c for different topological actions, with vortex suppressing 
parameter A = 0, 2 and 4, based on fits to the function 

This suggests a phase diagram as sketched in Figure [TJ We expect the endpoint 
of the transition line to be located at (A, 5) = (+oo, vr), see Section 4. 



7t/2 



massive phase 



massless phase 



Figure 1: A schematic illustration of the phase diagram, as expected based on the 
results for 5c{X) in TableUl and anticipating the outcome of Section 4- 
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3.2 Continuum Limit in the Massive Phase 



In order to investigate the continuum limit in the massive phase, we consider the 
step-2 SSF 3 

S(2, u, a/L) = ALm{2L) . (3.2) 

Here u = 2Lm{L), and m(L) is the size-dependent mass gap. Based on the exact 
S-matrix of the sine-Gordon Model, the SSF has been worked out analytically in 
the continuum limit a{2,u) 



Using the_ standard action. 
We mentioned 



S(2,M,a/L^0) [25|,l26| 
this analytic result has been confirmed in numerical simulations 16 
before that here even the cut-off effects were predicted to have universal features. 
This refers to a lattice SSF of the form 



S(2,M,a/L) = a{2,u) + 



[log(e/a) + U]' 



+ 



log^(e/a) 



(3.3) 



where ^ = l/m{L — )■ oo) is the correlation length in infinite volume. 

Figure |2] illustrates the cut-off effects of the SSF at m = 3.0038 for the stan- 
dard action, and for the constraint topological action with the vortex suppression 
parameter A = , 2 or 4. The curves are fits to eq. fl3.3l). w here we have inserted 
the analytically predicted continuum SSF a{2,u) of Ref. 16|. As in the case of the 
2d 0(3) Model, the standard action approaches the continuum limit from above, 
whereas the topological actions approach it from below. The continuum limit of the 
step-2 SSF at u = 3.0038 amounts to a{2, u) = 4.3895, and the cut-off parameter 
c = 2.618 [16] was supposed to be universal. In Table |2]we list our results, obtained 
by fitting the parameters a{2,u), c and U to the lattice data. The data for the 
standard action are taken from Ref. p^], where only cr(2, u) and c were fitted, since 
f/ = 1.3 is known from perturbation theory. These results indicate that all different 
actions converge to this continuum limit. For the topological actions, however, the 
fits yield negative and A-dependent values for c. This suggests that both parameters, 
U and c, depend on the lattice action and are therefore not universal. 





a{2,u) 


c 


U 


xVd.o.f 


standard action 


4.40(2) 


2.4(6) 


1.3 


0.84 


A = 


4.421(28) 


-4.0(3.6) 


4.1(2.2) 


0.15 


A = 2 


4.427(23) 


-5.26(45) 


-0.31(8) 


2.51 


A = 4 


4.71(25) 


-21(9) 


-0.87(60) 


0.23 



Table 2: Fitting results for the cut-off effects of the SSF in eq. Ii3.3\) for various 
lattice actions. The data for the standard action are taken from Ref. tldj : they 
were obtained by fitting a{2, u) and c, whereas U is known perturbatively. For the 
topological actions at X = 0, 2 and 4, we fitted cr(2,u), c and U . 

To illustrate the compatibility of our data with the analytic prediction, we follow 
Ref. [lil and plot in Figure[3]the same data as a function of (t/ + log(^/a))~^, where 
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Figure 2: Cut-off effects of the SSF S(2,u,a/L) at u = 3.0038 for the standard 
action (data from Ref. fid]), for the topological action without vortex suppression, 
A = 0, and with explicit vortex suppression, for A = 2 and A = 4. All curves are fits 
to eq. ^3. 3\) . where we insert the continuum limit ct(2,m) = 4.3895. 



^ is still the infinite volume correlation length, and U is a fitting parameter that 
differs for each action. In this plot we have again constrained the continuum limit 
of the SSF a{2,u) to its analytic prediction. 

3.3 Critical Behaviour in the Massless Phase 

In contrast to second order phase transitions, only two critical exponents — com- 
monly denoted as rj and S — are defined in the conventional way also for the essential 
phase transition, which occurs in this model, cf. eq. (13. ip . Based on Renormalisation 
Group techniques, their values have been predicted to coincide with the correspond- 



ing exponents in the 2d Ising Model [27|]. Here we focus on the exponent rj, and its 
property to characterise the divergence of the magnetic susceptibility x- 
The corresponding relation and the predicted critical value of rj are 

^\^\ f C"^"^ massive phase _ ,, 




^ massless phase 

in a square volume V = L"^ . We now focus on the massless phase and insert the 
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continuum limit 
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standard action i- 
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---Q ' 
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0.02 0.04 0.06 0.08 0.1 0.12 0.14 

(U + log(^/a))"2 



Figure 3: Cut-off effects of the step-2 SSF E(2, u,a/L) atu = 3.0038 for the standard 
action (data from Ref. [16]), as well as for the topological action without vortex 
suppression, A = 0, and with explicit vortex suppression, for A = 2 and A = 4. The 
curves are fits to eq. i\3. 3\} . where we have inserted the continuum limit a{2,u) = 
4.3895. The values on the horizontal axis depend on the fitting parameter U , which 
is different for each action. Note that the plots in Figures\^ and\^ contain (mostly 
invisible) error bars in both directions. 



measured values of x iiito the formula 

where C is the proportionality constant of eq. (13. 4p . At least within the massless 
phase, i.e. for 5 < 5c, it should be possible to find a constant C, which makes the 



results for rj in different volumes coincide to a good approximation [31 . 

Figure H] shows our results for A = 0, 2, 4, and L = 128, . . . , 1024, with the 
optimal choice for the constant C at each A. We see that the qualitative prediction 
of a coincidence of the rj values in different volumes, up to some limiting angle 
(5iirnit , is well confirmed. One is now tempted to interpret ^umit as an estimate for 



Sc [31(1 . Table [3] shows that these values match the expected magnitude, but they 
are significantly higher than the precise results for 6c , given in Table [H Hence the 
coincidence of r] persists even in some (narrow) region of the massive phase (although 
eq. (13. 5p does not apply anymore). 
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A 





2 


4 


"^limit 


1.825(5) 


1.93(1) 


2.17(5) 


ri(6r.) based on eq. (13.51) 


0.255(2) 


0.278(2) 


0.301(1) 



Table 3: Results for the limiting angle ^umit for the coincidence of the rj values in 
different volumes, and for the critical exponent rj^ obtained from relation ^3.5\) . 



If we naively extract the r^-values at 5c, we obtain results for rjc, which are again 
in the predicted magnitude, but without a satisfactory precision, see Table [31 The 
rjc values determined by this simple method tend to be too large, in particular for 
sizable A values. 

Similar problems are notorious in numerical studies of the standard action, the 
Villain action, the Step Model etcJl The situation improves as one includes a loga- 
rithmic correction to the finite size behaviour of x, which has also been elaborated 
analytically in Ref. (27l|FI 

XocL2-''(lnL)-2'- , re = -1/16. (3.6) 

Much of the literature that dealt with conventional lattice actions focused on at- 



tempts to evaluate the critical exponent j22|, |32|-|38|| . Its numerical measurement 



is extremely difficult, as expected for a small exponent of a logarithmic term. An 



overview of the results on this long-standing issue is given in Ref. [39[. Only in 
2005 Hasenbusch reported a value which seems to confirm the prediction decently, 
Tc = —0.056(7) [ioj. However, in his study of the standard action on lattices up to 
size L = 2048, Hasenbusch had to fix = 1/4 as an input, and to introduce yet 
another free parameter by extending the logarithmic factor to (const. + InL)"^*^. 

We first try to estimate the exponents rj^ and by fitting our data on lattice sizes 
L = 128, . . . , 1024 measured at 6 angles slightly above and below 6c ■ The fits have 
a good quality, and the results are given in Table HI The theoretical value r^c = 1/4 
is reproduced well at A = and approximately at A = 2. However, at A = 4 we 
obtain an r]c value which is clearly too large. Nevertheless this is compatible with 
the scenario that the topological actions considered here are in the BKT universality 
class, and that the finite size effects are amplified for increasing A — in qualitative 
agreement with the observations of Subsection l3.2[ Since a sizable A value suppresses 
the vortex density, it takes a very large volume to provide a sufficient number of 
vortices to drive an (approximate) BKT transition — in line with the picture of 



Ref. 14 



Regarding the logarithmic term in eq. (13. 6p . we do obtain small exponents of 
\rc\ = C(O.l) or below, but within this magnitude we cannot reproduce of the exact 
prediction. 



direct consideration of the correlation function {cxf^x+r) r^^ is plagued with even worse 
practical problems. 

•^On the other hand, this logarithmic correction term hardly affects the plots in Figure S) 
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Figure 4: The dependence of the exponent t], according to eq. i\3. 5\} . on the constraint 
angle 6 at X = 0, 2 and 4. 
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A 


5 


n 

1 


r 







1.76 
1.78 


0.2563(66) 
0.2446(68) 


-0.016(19) 
-0.034(19) 


0.022 
0.05 


2 


1.86 
1.87 


0.255(26) 
0.2558(15) 


0.060(74) 
0.070(14) 


0.47 
0.11 


4 


1.92 
1.94 


0.366(11) 
0.317(25) 


-0.194(32) 
-0.0470(66) 


0.087 
0.012 



Table 4: Results for the determination of the exponents r] and r in eq. liS. 6\) . by 
fitting our data at L = 128, 256, 512 and 1024, in the vicinity of the critical points. 



Motivated by the strong finite size effects in this model, Refs. 22|, |3J] worked 



out even a sub-leading logarithmic correction, which extends ansatz (13. 6|) to 

X = L-^(lnL)--(a, + a,i^), (3.7) 

where ai and 02 are constants. We add fitting results to this extended formula, based 
on our data measured at 6c with fixed exponents rjc = 1/4, Tc = —1/16, such that 
only ai, 02 are free parameters. Table |5] and Figure [5] show that these data match 
this form accurately for A = 0, 2 and 4, if we consider some range with L > 1280 
This observation provides satisfactory evidence that the behaviour in these points is 
compatible with the BKT characteristics, so that the topological actions do belong 
to the standard universality class, in agreement with Subsection 13.21 We assume 
this behaviour to persist for all points on the transition line with < A < 4. The 
limit A — i- +00 will be addressed in the next section. 



4 Continuum Limit of the pure Vortex Suppres- 
sion Action 

We now investigate the vortex suppression action without an angle constraint (which 
corresponds to 5 = vr). Thus we consider the upper axis in the phase diagram of 
Figure [TJ We have determined the infinite volume correlation length ^ as a function 
of the vortex suppression parameter A on lattice sizes up to = 2000 x 2000. The 
results can be fitted well to the function 

e(A) = aexp(6A) , (4.1) 

where a and b are fitting parameters, see Figure |6l This suggests that the critical 
value is at A = +00, as we anticipated in Figure [TJ This limit can be viewed as a 
plaquette constraint action. 

^ Table [S] also shows that the ratio |a2/ai| increases rapidly with A. Hence the supposedly 
sub-leading term in eq. (j3.7l) dominates more and more, which is consistent with the previous 
observation that finite size effects are very strong at A = 4. 
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512 


0.1618(46) 


1.488(16) 


0.002 






32 


-0.073(18) 


1.999(52) 


25.45 






64 


-0.0232(68) 


1.847(20) 


1.548 


A = 4, 


(5 = 1.9361 


128 


-0.0005(36) 


1.774(12) 


0.131 






256 


0.0128(53) 


1.730(18) 


0.050 






512 


0.0406(98) 


1.634(33) 


0.014 



Table 5: Fitting results for the data at the critical angle 6c, in the range L^i^ to 
LmsLx = 4096. We fit the magnetic susceptibility x to eq. ( [5*. 7| ), with the predicted 
critical exponents rjc = 1/4, Vc = —1/16. For Lmin > 128 the fits work very well, 
which confirms the compatibility of our data with the critical behaviour of the BKT 
universality class of the 2d XY Model. 



Studying the transition by measuring the step-2 SSF 

a(2,u) = hm S(2,M,a/L) (4.2) 

a->-0 

(cf. Subsection l3.2l) confronts us with an additional limitation. The numerical results 
show that for this action the finite size effects constrain the finite volume correlation 
length to ^{L) < 0.4 L. This restricts the range of the variable u = 2m{L)L = 
2L/^{L) to a regime u > 5.0. 

A restriction of this kind is natural in models with discrete energy eigenvalues 
oc 1/L in a UV conformal limit j4l|Jf| Also for the standard action in the 2d XY 
Model there is an upper bound 

'-^^'.^O(^) (4.3) 



TT \logL 



in the massive phase, see e.g. Ref. [40| and references therein. Qualitatively, such an 
upper bound can be understood using inequalities for ferromagnetic systems. This 
is briefly discussed in Appendix B. 

We can still measure the step-2 SSF for u sufficiently large, for instance u = 
2m{L)L = 6, and try to fit the cut-off behaviour with the function from eq. (13. 3p . 



^This is the ordinary case; asymptotically free theories (in the usual sense) are the exception, 
where any u G is possible. 



13 



x = o 



le+06 7 



100000 7 



10000 7 



1000 




100 



le+06 It 
100000 
10000 
1000 



100 



le+06 



100000 7 



1000 



L 

A, = 2 







' ' 1 


r'' 












,e'' 

3''' 










. . 1 





1000 



L 

A, = 4 




Figure 5: Numerical results for the susceptibility x> rneasured at the critical points 
for A = 0, 2 and A, on lattices of size L = 128, . . . , 4096. The fits refer to eq. ( [5*. 7| j 
with fixed exponents r]c = 1/4, Vc = —1/16, and ai, a2 as free parameters. Here 
and in Table we see that these fits are accurate in all three cases, confirming the 
compatibility of our data with a BKT phase transition. 
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350 




Figure 6: Correlation length ^ on large lattices as a function of the vortex suppression 
parameter A. The fit to the exponential function (with a = 0.492, b = 0.729, 

X^/d.o.f. = 0.290y) indicates an essential phase transition at X = +oo. 



which describes the continuum hmit at a BKT point. This fit, shown in Figure 
[71 works quite well. However, its continuum extrapolation cr(2,M)fit = 9.474(12), 
given in Table is rather far from the analytic BKT value of a{2,u) = 11.5314 
|4l| (which is close to 2u = 12). This suggests that the endpoint of the transition 
line does not represent a BKT phase transition. Indeed, this point is specific in the 
sense that one cannot cross it (on the axis 6 = n). Moreover, this observation is 
fully consistent with the established pictures of vortices driving the BKT transition 



IJ], so it cannot occur in the absence of vortices. 



xVd.o.f. 


c 


u 


0"(2,M)fit 


a{2,u) 


0.72 


-0.11(10) 


0.44(68) 


9.474(12) 


11.5314 



Table 6: Fitting result for the cut-off effects of the SSF S(2,u,a/L) at u = 6, 
according to eq. hS. 5]) . for the pure vortex suppression action. The fitted continuum 
extrapolation cr(2,M)fit does not agree with the BKT value a{2,u). 

In the 3d XY Model, the analogous point (A = +oo, with no other restriction) 



42 . Also in that case the observation of the transverse 



has been studied in Ref. 
susceptibility oc L°-^ did not clarify the properties of this vortex-free case. 
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9.49 i-r 



9.48 -; 
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(U + log(5/a))-2 

Figure 7: Numerical data for the step-2 SSF E(2,u,a/L) at u = 6, for the pure 
vortex suppression action, fitted to eq. \3. 'J\) . The parameters are given in Tahle\^ 
The continuum extrapolation cr(2,-u)fit does not agree with the BKT value. 

5 Concluding Discussion 

In this paper, we have investigated topological lattice actions for the 2d XY Model. 
At the classical level, these actions do not define a proper field theory, and pertur- 
bation theory is not applicable. 

In order to efficiently simulate topological actions, we have employed variants of 
the Wolff cluster algorithm. Its application to the constraint action is straightfor- 
ward, and for the vortex suppression action a generalisation to 4-spin interactions 
has been developed and applied successfully, see Appendix A. 

Despite its classical deffciencies, just as in the 2d 0(3) Model, we found that — 
up to moderate vortex suppression — topological actions yield the correct quantum 
continuum limit, which is here associated with the BKT phase transition. This 
includes in particular topological actions where vortices do not cost any energy. 
This observation is remarkable in light of attempts to derive the critical line from 
the energy requirement for isolated vortices. 

Specifically, in the massive phase, just as for the standard lattice action, the 
continuum limit is related to the sine-Gordon Model. In the massless phase we have 
verified the usual BKT behaviour of the critical exponent rj^ — and of its logarithmic 
correction term — with the topological actions. Our study demonstrates again the 
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immense robustness of universality in quantum field theory, which does not rely on 
classical concepts. 

An exception is the endpoint of this critical line, which seems to be located at 
(A, (5) = (+00, vr). The extrapolation to this point — which represents a plaquette 
constraint action — does not coincide with the BKT behaviour. This agrees with 
the established picture that vortices (which are completely eliminated at this point) 



are required to arrange for a BKT transition [14 



For comparison, we mention the case of the so-called Extended XY Model, with 



the lattice action 43 



cos 



2g 



((<^x-<^y)/2) 



(5.1) 



in the notation of eq. (12. ip . and with g > 0. For g = 1 it is equivalent to the standard 
action, but increasing q leads to a more and more narrow potential well for ip^ — (py, 
with width ^ tt / ^/q. The motivation was also an explicit vortex suppression; in the 
quadratic approximation to the potential they cost energy ^ 

Due to the gradual suppression of the vortices for increasing exponents (even 
without fully excluding them), Ref. 43| predicted the phase transition to turn into 
first order above some value of q, so it would match the behaviour which is observed 
experimentally for melting films of noble gases adsorbed on graphite. 

This Extended XY Model has been investigated in numerous papers. The es- 
sential BKT phase transition is observed at low values of q, and for some time the 
conjectured first order transition at lar ge q was controversial. However, it is now 



well confirmed numerically at o> 8 j4J, |45i. Moreover, an analytical proof for this 



conjecture was given in Ref. j46|. Ref. j47| added a vortex eliminating term with 
A — )■ +00 also in this case. No phase transition was observed at finite /3, hence 
the authors concluded that not only the BKT transition, but also the first order 
transition at large q is driven by vortices. 

In contrast, for the Step Model no non-BKT phase transition has ever been found, 
and for the topological lattice actions we do not observe any finite order transition in 
the 5-X phase diagram either. However, a change to first order along the transition 
line — at some large value of A — is conceivable in our case as well (that would not 
contradict universality). If this occurs as in the Extended XY Model, then it should 
change again at the endpoint, according to Ref. [47 . 



In any case, the characteristics of the transition at the vortex-free endpoint is an 
open question, to be explored in the future. 
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A Cluster Algorithm for the Vortex Suppression 
Action 

The algorithm for the vortex suppression angle constraint action is based on the 
Wolff cluster algorithm In the single cluster variant, each cluster update begins 
with the selection of an initial spin as a seed for cluster growth, and with the 
choice of a reflection line (a reflection hyper- plane in general 0{N) Models), which 
is perpendicular to the randomly selected unit vector r. Starting with this seed, 
some spins e*c are combined to a cluster, which are then collectively reflected — or 
flipped — to the new spin orientations Cx' = Cx — 2(f ■ e^) f. Spins may be put in 
the same cluster due to the nearest-neighbour angle constraint, or due to the vortex 
suppression plaquette interaction. Two nearest-neighbour spins Cx and Cy are always 
put in the same cluster if the flip of Cx to Cx ' (without flipping Cy) would lead to a 
relative angle between Cx ' and Sy beyond the constraint angle 5. 

The cluster rules implied by the vortex suppression four-spin plaquette action 
are more complicated. Let us consider the spins Cx^ at the four corners xi, X2, 
and X4 of a plaquette □, as well as their reflection partners e^;. '. Depending on 
whether a spin is flipped or not, there are 16 possible spin configurations on the 
given plaquette. Each one has a Boltzmann weight exp(— Altinl), depending on the 
vortex number \vu\ of the corresponding spin configuration. Since the reflection of 
all four spins on a plaquette □ just changes the sign of the vortex charge vu, each of 
the 16 spin configurations has a total reflection partner with the same Boltzmann 
weight. We can thus limit the discussion to 8 pairs of configurations. We distinguish 
two qualitatively different cases: 

1. In this simple case the vortex number is always zero, irrespective of whether 
any spin is flipped or not. Hence all 16 spin configurations have the same 
Boltzmann weight 1. Based on the vortex suppression action, there is no need 
to put any of these four spins in a common cluster. 

2. The second case can be characterised as follows: when all spins are flipped to 
the same side of the reflection line, the vortex number is necessarily zero. We 
denote this spin configuration as the "reference configuration" . When each of 
the spins is individually flipped (without flipping any other spins), there are 
two spins whose flip generates a vortex (or an anti- vortex). We denote these 
two as the "active spins" . It turns out that the simultaneous flip of two spins 
(starting out of the reference configuration) generates a vortex only if exactly 
one of the two spins is active. If both or none of the two flipped spins are active, 
no vortex is generated. If three or four spins are flipped simultaneously, one 
just generates the total reflection partners of the previously discussed cases. 
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This gives rise to the following cluster formation rule. If the two active spins 
are on the same side of the reflection line, they are put into the same clus- 
ter with probability 1 — exp(— A), otherwise they remain independent. The 
other spins are not affected by the vortex suppression action on this plaquette 
and remain independent. Still, preliminarily independent spins may finally 
become members of the cluster due to the angle constraint, or due to the vor- 
tex suppression action on a neighbouring plaquette]^ For the efficiency of the 
algorithm it is essential that spins are put in the same cluster only if they 
are on the same side of the reflection line. This prevents the clusters from 
becoming unphysically large (their linear size should be of 0{C,)). 

This algorithm obeys detailed balance. In particular, when a plaquette carries 
a vortex (and thus has the Boltzmann weight exp(— A|fn|) = exp(— A)), the two 
active spins are not put in the same cluster (with probability w = 1), because 
they are then necessarily on two different sides of the reflection line. On the 
other hand, if the two active spins are on the same side of the reflection line, 
= and the Boltzmann weight is exp(— A|fg|) = 1. In that case, the two 
active spins are put in the same cluster with probability l — w' = 1 — exp(— A), 
while they remain independent with probability w' = exp(— A). Only in the 
latter case, the two active spins may not belong to the same cluster, and are 
thus flipped independently, which again results in the creation of a vortex. 
Hence the detailed balance relation connecting the two configurations reads 

exp(— A|fn|) w = exp(— A) = exp(— A|fg|) w' . (A.l) 

As we have explicitly verified in an extensive computer search, other cases do not 
exist. Once spins have been put together in the same cluster (due to the nearest- 
neighbour angle constraint action, and/or due to the vortex suppression plaquette 
interaction), all spins in the cluster are simultaneously flipped to ' ■ Then a new 
random site is selected as a seed for cluster growth, along with a new unit vector r, 
and the entire procedure is repeated. 

As an alternative to this single-cluster algorithm, we also employed a multi- 
cluster algorithm, which constructs all clusters in a spin configuration and flips each 
of them with a probability of 1/2. Then the subtleties explained in footnote E] do 
not occur. 

An additional virtue of cluster algorithms is the applicability of improved esti- 
mators. For the variant that updates the vortex suppression angle constraint action, 
the improved estimators — for example for the correlation function and the suscep- 
tibility — work exactly as in the original Wolff algorithm j2| . 

^It should be noted that two active spins that are tied together in the same cluster may actually 
end up not to belong to the single cluster that is currently being built. In any case, one must keep 
track of the plaquettes on which a decision based on \vu\ has already been taken, and one must 
stick to that decision when this plaquette is visited again, in the process of identifying the cluster 
members. 



19 



B On Inequalities for Ferromagnetic Systems 



Consider the standard action (12.11) on a long strip with N = L/a sites on a time- 
shce. (We take below a = 1 for simplicity.) Making the ferromagnetic coupling 
anisotropic, (3 — )■ {(3x,/3t), we increase (3x oo while keeping (3t = [3 constant. This 
way the 2d system turns into a Id chain with /3' = N [5. By increasing a /3-parameter 
in a ferromagnetic system one might expect that the correlation length could only 
grow. This implies a lower bound for the correlation length in the original model 
(with isotropic coupling) 

m N) < UNP) = 2f3N + 0{1) , (B.l) 
where G(/3') is the correlation length for the Id chain. 



According to Ginibre's Theorem [48| this intuitive argument indeed holds for the 



standard action, and for a large class of further actions specified in Ref. [48 . 



Surprisingly, for slightly more complicated actions this inequality does not hold. 
Consider the nearest neighbour action with the action density 

s(e, e") = /3(1 - e ■ e') + 7(1 - e ■ e'f + Sconstr(e ■ e' - cos 6) , (B.2) 

where the last term describes the constraint e-e' > cos 5. One can make the system 
"more ferromagnetic" by increasing or 73,, or by decreasing 6^^- Taking again the 
Id limit (say, by (3^ — > 00) one would naively expect 

aP,j,S;N)<UN(3,N^,S) . (B.3) 

This, however, cannot be true, since one has ^i(0, A^7, vr) oc y/W^ for N'-j — )■ 00 and 
^1(0,0,5) oc 1/(5^ for (5 — )■ 0, while the left-hand-side increases oc iV in the massless 
phase. (Of course, the action (1B.2p does not satisfy the conditions of Ginibre's 
Theorem.) 



/3x 


A 




It 




0.0 


0.0 


1.0 


1.0 


8.8316 


0.0 


0.0 


1.1 


1.0 


8.7980 


1.0 


1.0 


1.0 


1.0 


13.9148 


1.0 


1.0 


1.1 


1.0 


13.9008 


1.1 


1.0 


1.0 


1.0 


13.8877 


1.0 


1.0 


0.0 


0.0 


6.4263 


1.1 


1.0 


0.0 


0.0 


6.4553 



Table 7: The change of the correlation length by increasing the spatial parameters 
Px or 72,. on a strip with N = 2 sites. The last pair of data refers to the standard 
action where Ginibre's Theorem applies, so that the intuitive expectation holds. 

In Table [7] we illustrate this behaviour for the mixed action (without the con- 
straint, 6 = n), where the inequality is violated, and for the standard action (7 = 0, 
6 = it) where it holds. 
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